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SUMMARY 

Explicit solutions are obtained in parametric form for the magnetic field due to sur- 
face currents flowing on an infinitely long rectangular conductor. The well-known 
analogy between the two-dimensional magnetic potential problem and the corresponding 
electrostatic problem is exploited to obtain these solutions. The method employed is a 
previously developed Schwartz -Christoffel transformation that maps the complex space 
outside a rectangle into the complex space outside a unit circle. Explicit solutions for 
the magnetic field outside the rectangle are derived. 


INTRODUCTION 

The development of small electronic systems has stimulated considerable interest in 
the use of conductors with rectangular cross sections, such as thin films, tapes, etc. 
Application of these configurations requires a knowledge of the magnetic field created by 
the currents passing through the systems. In ordinary conductors where the current is 
uniformly distributed, the magnetic field can be obtained by straightforward application 
of the Biot-Savart law, but if the current only flows on the surface of the samples, such 
as with superconductors or with rapidly oscillating fields, minimum energy configura- 
tions must be considered, and the problem becomes more complex. In recent years, 
attempts to use a Green's function technique to solve for the current density in super- 
conductors (refs. 1 and 2) led to an integral equation that could not be integrated ex- 
plicitly. Solutions are generally obtained by numerical integration using successive ap- 
proximations. A simple alternative approach is used herein. This approach is based on 
the fact that a good approximation of the magnetic field configuration around the rectangle 
can be foimd at distances larger than the penetration or skin depth away from the con- 
ductor by assiiming that all currents flow on the surface of the conductor. With this as- 
sumption, a qualitative description of the current distribution can be obtained. This as- 
sumption is valid particularly for superconductors and normal conductors at high fre- 
quencies. 



For the case considered herein, infinitely long conductors of constant rectangular 
cross section with the current flowing only on the surface (e. g. , a perfect conductor), 
a close connection between some electrostatic and magnetic problems exists. Knowledge 
of the relation between the electric and magnetic field potential functions in these in- 
stances allows the use of techniques developed for solutions of Laplace's equation, such 
as complex transformations. Through an analogy between the corresponding potential 
functions, explicit equations are derived for the magnetic field and surface current dis- 
tribution by use of the conformal transformation of a rectangle into a unit circle (ref, 3), 
The resulting field is described by equations in which the parameters are the coordinates 
of an intermediate space. 


ANALYSIS 

Fundamental Mapping 

For a conductor of zero resistivity, no internal electric field can be maintained. 
Maxwell’s equation for a time-varying magnetic field V x E = - (0B/3t) then leads to the 
conclusion that the magnetic field S is a function of coordinates alone. In order to 
simulate the case of superconductors, the Meissner condition that B = 0 inside is used. 
This case also applies to normal conductors with high-frequency alternating currents. 

Since no magnetic or electric field is allowed to exist in the interior of this perfect con- 
ductor, the current density is also zero, as indicated by Maxwell's equation for a time- 
varying electric field V x B = J + e(0E/9t). At the boundary, however, surface currents 
and surface charges are possible. The tangential component of magnetic field is dis- 
continuous at the surface of the conductor by an amount proportional to the surface cur- 
rent per unit width. The normal component of magnetic field is continuous across the sur- 
face, and since it is zero inside the perfect conductor, it must go to zero outside the 
bovmdary. In analogy, the normal component of electric field is discontinuous at the sur- 
face by an amount proportional to the surface -charge density. If the perfect conductor is 
assumed to be sufficiently long so that the distribution of current is not affected by contacts, 
a two-dimensional analysis is possible. 

As shown in some texts on electromagnetic theory (ref. 4), an analogy exists be- 
tween the magnetic potential A due to a surface current and the electrostatic potential 
(p due to a surface -charge density. For a conductor viiose axis is taken in the direction 
of the z-axis, the basic solutions for each type of problem are as follows: 
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Electrostatic: 


Magnetic: 


(p =_L_ /* P(x*,yOdv» 


Mq 

n J If - 


A = kA, = k 

Z 47T 


( 1 ) 


(All symbols are defined in the appendix. ) Thus, if p_ and J_ are similar functions of 

the coordinates, then <p and A^ (as well as E and B) will have similar solutions pro- 
vided that <p and A^ satisfy similar boundary conditions. 

The necessary and sufficient conditions to define a unique potential <p in all space 

are 


^1 = ^2 constant 


and 


3n 


- 


3<jpj 

3n 


a 



where and ^2 the potential fxmctions in regions 1 and 2, respectively, /Cj and 
«2 are the dielectric constants, and a is the surface -charge density. 

The equivalent boundary conditions for A^ are derived as follows: If Bj and §2 
are the magnetic fields inside and outside the rectangle, respectively, and n is a unit 
vector perpendicular to the surface, the following equation results from V • B = 0: 


(B 2 - Bj) • n = 0 

where 



Then, since B = vA„ x k, 

2i 


(VAg 2 ^ ^ 2 X k) * n = 0 


which for the rectangular coordinate system leads to 


n X (VA^^ 2 ■ l) - 0 


and 


A„ n - A„ 1 = constant 

z, ^ Z, X 


The secondary boimdary condition on A„ is obtained as follows: 

z 


\ - 

nXl — 1 = J 


or 


Mg Mj, 


n X ( — VA, « X k - — VA^ 1 x£) = J 

Y2 ’ ^1 ’ / 


which for the rectangular coordinates leads to 




k = - J 


and finally 


1 ^^z,2 1 


lM2 


Ml dn 


k= - J 


Since both and <p obey Laplace's differential equation inside and outside and are 
subject to the same boundary conditions, the uniqueness theorem indicates that the solu- 
tions for the two potential fields are identical (ref. 5) hence, A^ = Cq<p where Cq is a 
constant used to obtain the proper units. Now 
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( 2 ) 


E = - grad (p = - — grad 
^0 


and 


B = CqIc X E 


(3) 


Equation (3) shows that B and E are perpendicular to each other at any point. Since 


3 ^ 

47t y I r - r* I ® J I ^ " ^ I 


dv’ 


(4) 


then 


Jz(x’,y’) 


^0 

— P(x'.y’) 


(5) 


where J_(x',y*) is a surface current density and p(x',y*) is a surface -charge density. 
Since the analysis is two-dimensional, is a current per unit width, and p(x',y’) is a 
charge per unit width per vmit length. Hence, two-dimensional ms^netic problems allow- 
ing surface currents only can be solved by using techniques common to two-dimensional 
electrostatic potential problems. One such method involves the use of a conformal trans- 
formation of complex spaces. This transformation facilitates the solution of Laplace's 
equation by mapping into a space of simpler geometry. 

In solving potential problems by use of the conformal transformation, it is often as- 
sumed that, once the transformation equations are derived, the solution is essentially 
complete. Certain parameters of the problem, such as the charge distribution at the sur- 
face and the shape of the equipotentials, may be found directly from the transformation 
equation; however, the dependence of the field and of the potential on distance and angular 
inclination from the source is still unknown. The procedure to follow then is to solve 
Laplace's equations for the simpler geometry and then map the solutions into the more 
complex geometry by using the Schwartz -Christoffel transformation equations. Specifi- 
cally, the actual field and potential for a desired location are obtained by mapping the 
proper points outside the simple conductor into the desired point and by solving for the 
associated field and potential. 
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Solutions for the Rectangular Conductor 


Consider an infinitely long conductor of uniform rectangular cross section with sides 
2a and 2b (fig. 1). Current flow through the system is assumed to exist only in an in- 
finitesimally thin layer at the surface. The potential aroimd the rectangle can be found 
with the use of the Schwartz -Christoff el transformation, A transformation equation was 
derived in reference 3 that maps the two-dimensional region outside a rectangle in the 
complex £ -space (fig. 1) into the two-dimensional region outside a unit circle in the com- 
plex t-space. The following differential equation is derived in reference 3 to perform 
the transformation: 


— = F’(2 cos 2t - 2m)^/^ (6) 

dr 


where t = e^^. The parametric solution of this equation x = f(r) and y = g(r) described 
in reference 3 is not in a form that is easily usable; however, since tables and computer 
programs of Legendre polynomials are now available, a more readily usable solution may 
be derived in the following manner. Integration of equation (6) leads to 


t = 2F’ 


-i(2n + l)cos 2t +. 2 sin 2t ip 


4 - (2n + 1)" 


2n + 1 


+ D 


(7) 


where Pjj(m) are the Legendre polynomials. The constant of integration D corresponds 
to a displacement of the circle to which the rectangle is mapped from the origin of the 
t-space, and for simplicity is set equal to zero. The interpretation of ix = cos 2 q! in 
the t-space is extremely important. Each point on the circumference of the circle given 
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by 1±Q! and Tr±a maps into a corner of the rectangle. Therefore, a is a constant 
whose value dictates the ratio of the sides of the rectangle, that is, b/a. The relation 
between [i and b/a is shown in reference 3 to be expressible in parametric form; 

b ^ E(k) - (1 - k^)K(k) 

^ E (Vl -k^ ) - k^K (Vl -k^) 


and 


11 = 1 


- 2k 


2 


This relation is plotted in figure 2. Solving equation (7) in terms of the components of 
T, with T = 6 - i/3, yields the following relation between the two complex planes: 


£ = 2F’ ^ P^(/i)e"(2n+l)^[cos(2n + 1)0 _ i sin(2n + 1)0] 


n=l 




sin 26 


4 - (2n + 1)' 


1^2 cosh 2/3 + (2n + l)sinh 2/3] 


i cos 29 
4 - (2n + 1)" 


1^2 sinh 20 + (2n + l)cosh 2/3] - — 


( 8 ) 
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where the complex space t can be broken into components: f = x + iy. The constant 
F', as derived in reference 3, is F' = ib/2L where L = E - k’K and K and E are 
elliptical integrals of the first and second kind, respectively. From equation (8) the 
parametric equations for x and y can be written as 


b 


X = — 


L 




'' 

4 cos(2n + 1)6 cos 20 


2 sinh 2^ + (2n + l)cosh 2/3 
4 - (2n + 1)^ 


+ sin(2n + 1)0 sin 20 


2 cosh 2/3 + (2n + l)sinh 2^ 

M cos(2n + 1)0 

4 - (2n + 1)^ 

2n+ 1 

J 


(9a) 


oo 


n=0 


- sin(2n + 1)0 cos 20 
V. 


2 sinh 2j3 + (2n + l)cosh 2/3 
4 - (2n + 1)^ 


+ cos(2n + 1)0 sin 20 


2 cosh 2/3 + (2n + l)sinh 2/3 

jLt sin(2n + 1)0 j 

4 - (2n + 1)^ 

2n + 1 J 


(9b) 


On the surface, the series expansion involving the Legendre polynomials cannot be used, 
and the differential transformation equation of reference 3 (that is, eq. (6)) must be 
integrated as follows: 


^ = F' f(2 cos 2 t - 2ju)^/^ , 
dr L 

F'(2 cos 20 - 2 m)^/^ 
d0 


( 10 ) 


( 11 ) 


Integration of this differential equation using elliptic forms leads to the following: 
For sides BC and AD: 


f = 2F' [e( 0, sin O') - cos^a F(0, sin 0 !)J 


(12a) 


where 


-a ^ 9 ^ a and n - a^d^v + a 
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(c) Square. 

Figures. - Concluded 


For sides AB and CD: 


t = -2iF' [e(0, cos a) - sln^a F(0, cos a)] 


(12b) 


where 


(X < 6 < TT - a and v + a ^ 6 ^ ~ a 

In these two equations, the constant of integration has been set equal to zero so as to 
keep the rectangle centered at the origin of the f-space. E(0,k) is the incomplete ellip- 
tic integral of the second kind, and F(0, k) is the incomplete elliptic integral of the first 
kind. Both can be found in the appendix of reference 6 (note that k = sin a or cos a). 
The fundamental net, which corresponds to the lines of constant magnetic potential and 
the lines of magnetic flux, can be plotted from equations (9a) and (9b). The lines of con- 
stant magnetic potential correspond to the family of constant and the lines of flux to 
the family of constant 6. The lines of constant magnetic potential for three rectangles 
(length -to -width ratios of 10:1, 5:1, and 1:1) are plotted in figure 3 around the corner of 
each rectangle for the same values of 
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Current Distribution 


The use of equation (5) and the charge distribution of reference 3 shows that the 
normalized distribution of surface current per imit width in the region - a < 0 < a, 
which maps into the side BC (fig. 1), is given by 



( 13 ) 


where 


4(a + b) 

In the region a < 0 :< n- - a that represents side AB in figure 1, the current per unit 
width becomes 



+ A V2 

7T \b / 

/ 2 2 
ysin 0 - sin a 


(14) 


The other two sides of the rectangle (AD, DC) have the same corresponding current dis- 
tribution because of the symmetry. The current per unit width becomes infinite at the 
corners (0 = ±a, n±a) because the magnetic field goes to infinity at the corners. The 
total current, however, remains finite. The profile of this current density distribution 
across the conductor shows qualitatively the behavior of real currents flowing through 
superconductors and through normal conductors at high frequencies. The important point 
is that the highest concentration of current is around the corners. 

This behavior of current per unit width at the corners is inherent in problems con- 
cerned with pure surface currents. In actual superconductors, where the current at the 
corners does not go to infinity, a deep penetration of the magnetic field into the sample 
occurs that causes the current to flow in a thicker layer along each corner. 
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Magnetic Field Distribution 


The magnetic field can be obtained explicitly from the relation between the complex 
functions involved. A complex potential function P(0, /3) is sought in the t -space and 
then transformed into the £-space. When P(0, |3) is obtained for the f -space, the mag- 
netic field is then obtained by the following transformation: 

i=-CQkXT (15) 

where T is a vector in which the real and imaginary parts of (dP/df)* are the com- 
ponents in the x and y directions, respectively: 



(16a) 


Ty = imag 



*' 


The complex potential ftmction is written as follows: 

P(0,^) =U(0,/3) +iV(0,/3) 


(16b) 


(17) 


where in the t-space, U(0,/3) is the magnetic potential, and V{9,0) represents the lines 
of magnetic flux. The complex solution of Laplace's equation in the t-space is 


P(0,/3) = In t = In 6 + ie (18) 

where t = 6 e^^ = e^^ = The potential solution for the unit circle of the t-space 

is U = In 6; V = 0 represents the stream function. The magnetic field due to the 
normalized current of equations (13) and (14) may now be fotmd from equations (6), (15), 
and (18): 


B = -G k X (— — — \ = G k X 
o \dt drdf/ ° 


LiF'(2 cos 2t - 2m)^^^J 


( 19 ) 


where 
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Mq2(e + b) 


Substitution of t = 0 - i/3 and fi = cos 2 a into equation (19) leads to 


B = G k X . _ 

o 1-/2 iF 


Letting 


[(cos 26 cosh 2/3 - cos 2o!) - i(sin 20 sinh 2/3)J^^^ 
[(cos 20 cosh 2^ - cos 2a)^ + (sin 20 sinh 2/3)^] ^ ^ 


( 20 ) 


and 


M = cos 20 cosh 2/3 - cos 2o! 


N = sin 20 sinh 2^ 


( 21 ) 


yields the following simplified expression for B: 


B = G^k X 


(M - iN) 


1/2 




-/2 iF 9 <, 1/2 

(M'^ + N^) J 


( 22 ) 


The f -space must be divided into two regions that depend on the value of M to deter- 
mine the value of (M - iN) . In region I, let M s: 0. Several algebraic manipulations 
then yield the following expressions for the magnetic field components: 


-i/2 G L 


x,I 


— [(cos 20 cosh 2/3 - cos 2o!)^ + (sin 20 sinh 2/3)^j 


-1/4 


X 


sinf^tan-l sin 20 sinh 2/3 ) 

\2 cos 20 cosh 2/3 - cos 2a J 


(23a) 
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B. 


y,i 


^ r 2 2l “1/^ 

|^(cos 26 cosh 2/3 - cos 2 q')‘^ + (sin 26 sinh 2/3)^J 


/ 1 . -1 sin 26 sinh 2/3 

cos 1 — tan 

\2 cos 26 cosh 2/3 - cos 2a j 


(23b) 


In region II let M ^ 0, The magnetic field components are obtained in region II in a 
similar way: 


®x,n = * b 


Go -i/2 L 

- [(cos 26 cosh 2/3 - cos 20!) + (sin 26 sinh 2j3)^J 


cos / — tan"^ 


sin 26 sinh 2/3 


-cos 26 cosh 2/3 + cos 2o!; 


(24a) 


G, V2 L 


o ^ r 2 91 

By jj = ^ [(cos 26 cosh 2/3 - cos 2o!) + (sin 26 sinh 2jB)^j 


21-1/4 


X 


/l , -1 sin 26 sinh 2j3 

sin - tan 

\2 -cos 26 cosh 2/3 + cos 20! 


(24b) 


The absolute value of the magnetic field is 


|b|2 = 


2g5l^ 


^ — [(cos 20 cosh 2^ - cos 2o!)^ + (sin 20 sinh 2/3)^] 


- 1/2 


The magnetic field is tangent to the equipotential lines at all points. Direct substitution 
into equations (23) and (24) gives for the surface of the conductor (/3 = 0) in region I, 

B^ I = 0 and By j J ; for any /3 / 0 at 0 = 0 and at 0 = tt, j = 0 and By j = |b | . 

In region II for /3 = 0, jj J and By jj = 0; for any /3 / 0 at’ 0 = tt/2 and at 

0 = Zit/2, B^ jj = |b I and By jj = 0. At the border between regions^ and n, 

(cos 20 cosh 2/3 - cos 2 q! = 0), B^ j = B^ jj and By j = By jj. At infinity the denominator 

dominates the expression, and both B^ and By go to zero. 

Due to the complexity of equations (9a) and (9b), the magnetic field B cannot be ex- 
pressed in terms of x and y, the real variables of the original space. When knowledge 
of the magnetic field strength is sought for a given region in space, two ways of obtaining 
an answer are available. First, a table of corresponding values between x, y, and 0, jS 
may be calculated from equations (9a) and (9b) for the particular width -to -length ratio in 
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(b) y-coordinate against p. 

Figured - Coordinates x and y against p for rectangle of length -to-width 
ratio of 10. (0j and 02 are shown for a point (5, 9). ) 

question, and § can be found from these values of 6 and j8 from equations (23) and (24). 
The second method is more easily accomplished and is demonstrated as follows. Due to 
the symmetry of the rectangle, only the first quadrant 0 ^ 6 ir/2 need be considered. 
Points outside the first quadrant can easily be rotated into the first quadrant by sub- 
tracting nir/2, with n depending on the location of the point. Therefore, plots of x 
and y against ^ for a family of constant 0's between zero and tt/2 are obtained from 
equations (9a) and (9b) (fig. 4). For a given point, R(x, y), two values of 6 (e.g. , 0^ 
and 02) are picked from the x and y graphs such that the corresponding /3 values to 
each 0j and 02 are nearly equal. This is shown in figure 4 for the point R(5, 9). The 
values obtained for 0^ and 02 are 1. 04 and 1. 17, respectively. Once 0j^ and 02 are 
obtained, new plots of x and y against /3 can be calculated for a family of constant 
0's between 0^ and 02. New limits on 0 are obtained, and the process is repeated 
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until complex point R'(®) gives a point R'(x', y*) within the required accuracy of the 
original point R(x, y). Usually, the process need not be repeated many times, since, in 
practical applications, values of magnetic field are only sought for given areas, such as 
locations of wires or circuit components, rather than mathematical points. 


CONCLUDING REAAARKS 

The magnetic field of a rectangular conductor, whose current flows on the surface, 
was obtained by the use of a complex transformation and the corresponding two- 
dimensional electrostatic potential solution. The results can be applied to any rectangular 
conductor whose cross-sectional dimensions are large compared with the magnetic field 
penetration distance so that the induced current can be assruned to flow on the surface. 
Such systems are waveguides, good conductors with alternating currents, and super- 
conductors whose cross-sectional dimensions are much larger than the London penetra- 
tion depth. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 24, 1966, 

129-02-05-02-22. 
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APPENDFX - SYMBOLS 


A 

magnetic vector potential 

Jz 

component of surface current 

Az 

component of A along z - 
coordinate axis 


density along z -coordinate 
axis 

a 

B 

width of rectangular conductor 
magnetic field induction 

K(k) 

complete elliptic integral of 
first kind 

®x,I 

component of B along x- 

K» 

K(k’) 


coordinate axis in region I 

k 

modulus of Jacobian elliptic 
functions and integrals, 

®x,n 

component of B along x- 
coordinate axis in region n 


sin a 

B , 

comi)onent of B along y- 
coordinate axis of region I 

k’ 

complementary modulus, 

y,i 


Vl -k^ 

®y.n 

component of B aloi^ y- 
coordinate axis in region n 

£ 

imit vector along z -coordinate 
axis 

b 

thickness of rectangular con- 

L 

E(k) - k'^K 


ductor 

L' 

E(k’) - k^K’ 

^0 

1 Wb/(V)(m) 

M 

cos 20 cosh 2/3 - cos 2 q! 

D 

constant of integration 

N 

sin 20 sinh 2/3 

dv’ 

volume element of conductor 

n 

imit vector perpendicular to 

E 

electrostatic field 


surface 

E(k) 

complete elliptic integral of 2nd 

p(0,^) 

complex potential function 


kind 

Pn(M) 

Legendre polynomial 

E' 

F’ 

E(k’) 

1 _ ib 

- — C CSC or = — 

R(x,y) 

sample point in calculation of 
transformation 


4 2L 

r 

radial coordinate of field 

Go 

Mq 2(a + b) 


pomt, r = X + y + z 


If 

r' 

radial coordinate of source 
point, r*^ = x'^ + y'^ + z'^ 

I 

total current passing through 

f 


conductor 

vector whose components are 
real and imaginary parts of 

J 

surface current density in the 
conductor 


(dP/dS)* 
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t complex space where con- 

ductor is circular, 

t = e^^ = 6e^^ 

X, y, z coordinate axes of t -space 

a angle in t -space correspond- 

ing to corners of rectangle 

/3 In 6 

6 radial vector in t -space 

e permitivity 

permitivity of free space 

^ complex space where conduc- 

tor is rectangular, 
t = X + iy 

6 angle in t -space 

K|,K 2 dielectric constants 

IX cos 2a 

magnetic permeability of free 
space 


ix^, 1 X 2 magnetic permeability for 
regions 1 and 2 

p(x’,y’) electrostatic charge distri 
bution 

O’ surface -charge density 

T In t = 0 - i^ 

cp electrostatic potential 

Subscripts: 

1, 2 regions on either side of 

boundary between two 
materials 

I region I, 

cos 26 cosh 2)3 

- cos 2 q! ^ 0 

II region II, 

cos 20 cosh 2/3 

- cos 2o! < 0 

Superscript: 

* denotes complex conjugate 

functions 
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